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^ ■ Abstract 

O . An accurate simulation of Green's function and self-energy function of non-interacting electrons 

in disordered graphenes are performed. Fundamental physical quantities such as the elastic relax- 

K* ' ation time Te, the phase velocity Vp, and the group velocity Vg are evaluated. New features around 
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Introduction. — Graphene, a single layer of graphite, has been intensively studied in recent 
years because of many intriguing transport properties [l|, |2j. Examples include minimum 
conductivity and linear carrier density dependence of conductivity 3|, |J] . Electrons in an 
ideal graphene are governed by the relativistic massless Dirac equation and exhibit a linear 
dispersion relation in the vicinity of the Dirac point and zero density of state at the Dirac 
point p, |6| . Among many relativistic effects, the Klein paradox 7\ is arguably one of the 
most important effects that differentiates Dirac electrons [8] from the Schrodinger electrons 
in disordered systems. 

A deep understanding of Dirac electrons in disordered graphenes requires extended knowl- 
edge on the self-energy function in order to extract such fundamental physical quantities as 
the phase velocity Vp, the group velocity Vg, and the elastic relaxation time t^. However, it is 
known |9| that accurate and reliable calculations are quite difficult and nontrivial. For this 
reason, various approximations have been employed in different theoretical studies. So far, 
almost all calculations 10, lll(| concerning electron properties in disordered graphenes were 
performed without fully considering disorder effects. The wave-nature of Dirac electrons is 
more pronounced near the Dirac point because of very large electron wavelength there. It is 
known that interference due to multi-scattering leads to weak localization and the Ander- 
son localization in conventional disordered electron systems 12[|. Quantum interference also 
plays important roles in coherent wave propagation through quasi-random 13| and random 



media 
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15| . As it will be shown below, interference- induced the hybridization of Bloch 



states is essential in understanding the diffusion properties near the Dirac point in realistic 
disordered graphenes, although previous calculations [10|, lll| only captured the essential 
physics far away from the Dirac point. 

In this Letter, we present a systematic method to exactly calculate Green's function and 
the self-energy of large size disordered graphenes. We extract accurate Xg, Vp, and Vg values 
from the spectral function A(k, E) derived from the self-energy function. In comparison 
with the results from the self-consistent Born approximation (SCBA), it is found that Tg 
is overestimated by the SCBA in the strong disorder regime. We show that both Vp and 
Vg deviate significantly from their unrenormalized values and exhibit substantial energy 
dependence, especially near the Dirac point. The effective group velocity is larger than 
the effective phase velocity, but significantly lower than its unrenormalized value when the 
mixing of different Bloch states is dominant. Moreover, we generalize the Einstein relation 



to calculate the conductivity of the disordered graphene. 

Model- and-method. — 7r-electrons of undoped graphene can be modeled by a tight-binding 
Hamiltonian on a honeycomb lattice of two sites per unit cell, Hq = t ^ \i >< j\ + h.c, 
where t = —2.7 eV is the hopping energy. The corresponding eigenvalues and eigenstates 



of Hq near the Dirac point are [5|, ll7(|, respectively, -Efc,± = ^hvpk and |k± >= (|ky4 > 
_l_gj(^{k) |]^^ >)/2, where Vp is the unrenormalized Fermi velocity, and h is the Plank constant. 
A and B stand for A- and i?-sublattices. 0(k) is the polar angle of the momentum k, 
and |k74(5) >= .^ J2 e*'''''*(^'|ryi(B) >, where rA{B) is the position vector of the 
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(S) 



rA(,B) 



y4(i?)-lattice and Na{b) is the total A(i?)-lattice points. The plus (minus) sign denotes the 
conduction (valence) band. The Green function of clean graphene is, in a diagonal basis, 
G'o(k, E) = ^^,o+^_,^oj k+ X k+ | + g^,o+^^^^^J k- >< k- |. A weak point-like disorder is 
introduced through V = X]^«K >< "^1 with riimp randomly distributed impurity sites where 
the on-site energy ej of each impurity can take — Vq or Vq (measured in the unit of t) with 
equal probability. A dimensionless parameter a = "2^^^" i^" '^^^ ^^ used to characterize 
disorder strength. Here, Ac is the area of the unit cell. Our calculation shows that the 



physical quantities such as self-energy are only determined by the parameter a 18|. 
The ensemble-averaged Green function is defined as G(k±, E) 



<k±|-F^^ — Fj — rrlki >, where the bar means the ensemble average. It can be caL 



\ E+irj-Ho-V 



culated by using the well developed Lanczos recursive method |19l. |20[ |. In order to obtain 
an accurate ensemble-averaged Green's function near the Dirac point with high energy 
resolution 



2l[ |. a large sample containing A^ = L^ x Lj^ ^ 6.0 millions carbon atoms 
(2400x2400) is used in our simulation, where Lx{Ly) is the number of atoms in x{y) 
directions. The periodic boundary condition is used in our simulation in order to reduce 
the finite size effect. Thus, the wave vectors are kx = nx4:7i/3aLx and ky = nyATc/y/SaLy, 
where n^^y) is an integer and a is the lattice constant. 

Self- energy function. — Fig. 1 shows the calculated real (a) and imaginary (b) parts of the 
self-energy function for riimp/N = 10% and Vq = 0.5 (open squares) and 2.0 (open circles), 
respectively. The self-energy function is defined as usual, where T,(k,E) = GQ^{'k,E) — 
G^^(k,E) [23||. In principle, the self-energy function depends on energy E and wave vector 
k. However, our simulation finds that the self-energy function is not sensitive to wave vector 
k, while the one-particle Green's function depends on both E and k. This is not surprising 



since the scatterer size is much smaller than the electron wavelength near the Dirac point 
so that the inhomogeneous structure of disordered graphene can be described well by an 
effective homogeneous medium. This is also why the self-energy function is assumed to be 



k-independent in many perturbative calculations. Our finding validates this assumption 23|. 
Our calculation should be compared with widely used results from the self-consistent 
Born approximation (SCBA) that predicts 
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^ -2E/a-i2T,, \E\^T, 

-2a{E + i-Ka\E\)ln\^\ - i7ra|E|, \E\ > Tq 

where Fq = EcC^^^" (Ec is the cut-off energy) As shown in Fig. 1, SCBA results agree well 
with our exact self-energy function for very weak disorder (Vo=0.5 and rj-.vr.p/A^=10%)J25|. 
When the disorder strength increases several times (Vo=2.0 or a~0.07), the perturbative 
results can not capture the main features, especially near the Dirac point. The discrepancy 
is obvious for JmS(O) as shown in Fig. 1(b). luiL^E) at the Dirac point is Fq ~10^^t 
{Ec = vSt [6|) from the SCBA, while our exact value is /mS(0) ~10~^t. Thus, the true 
broadening of states at or near the Dirac point is much greater than what has been predicted 
by SCBA. This huge discrepancy can be attributed to the mixture of Bloch states caused 
by the impurities. Furthermore, the level-repulsion effect pushes all energy level toward the 
Dirac point so that the density of states at the Dirac point increases more in the presence of 
impurities. Thus, the impurities make the imaginary part of the self-energy function (directly 
associated with the density of states) at the Dirac point bigger. When the wavelength 
becomes short and the quantum interference as well as the Bloch state mixing are less 
important, IiriL{E) are determined by the disorder scattering. The difference between our 
exact simulation and that of the SCBA is small as shown in Fig. 1. 

Spectral function. — The single-particle spectral function relates to Green's function 
through A(k±,E) = -/mG(k±, E)/7r. Fig. 2(a) is A(k+,E) for ky = 0, I/q = 1, 
nimp/N=lQ%, and various kx (curves from the left to the right in the figure) ranging from 
0.0 (t2^=0) to 0.098 (ra^=56) (in unit of a^^). In the absence of disorder, the spectral func- 
tion Ao(k±,i?) is a delta function, reflecting that the wave vector k is a good quantum 
number and has all its spectral weight precisely at the energy E = E^j-. In the presence of 
disorder, the translational symmetry is broken and the spectral function is broadened due 
to the disorder scattering effect. The widths of the spectral function are given by ImTi{E) 
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FIG. 1: (Color online) (a) Real part of self-energy as a function of energy. (b)Imaginary part of 
self-energy as a function of energy. The open squares and open circles represent our numerical 
calculations for the disorder concentration riimp/N = 10% and Vq = 0.5 and Vq = 2.0, respectively. 
Dashed lines represent the SCBA results for the same disorder. Energy is measured in the units 
oft. 



-. Therefore, the elastic 



that measures the elastic relaxation lifetime t^, where Tg = _2im T,(E) ■ 
scattering relaxation time is akin to ImT,{E), and t^ around the Dirac point is mainly deter- 
mined by the Bloch state mixing and the level repulsion effect 26| . Far away from the Dirac 
point, the lifetime is mainly attributed to the disorder scattering. As shown in the inset of 
Fig. 2, lifetime become shorter as the wave vector increases. Physically, this is because the 
density of states p increases linearly with energy and the disorder scatte ring effects become 
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larger, which is qualitatively consistent with the prediction of SCBA 

Effective hand velocity. — Dirac electron propagation velocities, including the group 
velocity Vg and phase velocity fp, are also greatly modified by the disorder effects. These 
quantities relate to the shape of the dispersion relation that is the roots oi E — E^lk) — 



ReT,{E) = 0[27|, |28j. One can also extract the dispersion relation from the peak of the 
spectral function A{]i.,E) for a given k. Fig 3 (a) shows our exact dissipation curve Ef.ff{k) 
that is linear for very weak disorder {a < 0.01). However, when the disorder strength 
becomes significantly large {a ~ 0.1), the Eeff{k) becomes concave near the Dirac point, 
indicating the reduction in both the group and phase velocities. 

Fig. 3(b) shows A;— dependence of the group and phase velocities obtained from Vg = 
dEf.ff{k)/dk and Vp = Ef.ff/k. For the very weak disorder a = 0.01 as shown in the figure, 
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FIG. 2: (Color online) Single-particle spectral function A{k+^ E) plotted as a function of energy E 
at several k-points (from left to right): k = 0.000, 0.014, 0.028, 0.042, 0.056, 0.070, 0.084, 0.098 
(or Tia; = ~ 56) along the kx direction. The model parameters are Uimp/N = 10% and Vq = 
1.0. Inset: The energy dependence of single-particle relaxation time Te for Vq = 0.5 ~ 2.0 and 

Uimp/N = 10%. 

Vg and Vp are not too much different from the unrenomahzed velocity v^p. Their values 
are reduced by 5% in comparison with v^p. When the disorder strength increase several 
times a = 0.07, the renormalized Vg is higher than Vp and both Vg and Vp are reduced by a 
large percentage near the Dirac point. This shows that disorder not only renormalizes the 
Dirac electron velocities, but also changes the linear dispersion relation. The fact of large 
reductions of velocities at the Dirac point indicates that Dirac electrons near the Dirac point 
are more sensitive to the disorder. 

One can use the renormalization factor Z defined as Z = {1 ^qe )~^ [29] to measure 

the effect of disorder on electronic structure. Its value equals the ratio Vg/vp. As shown 
in Fig. 3(b), Z is very close to 1.0 for very weak disorder. In this regime, the Bloch state 
is still a good starting point for understanding disordered graphene, while the transport 
properties are expected to be described by the quasi-classical Boltzmann theory. When the 
disorder strength increases several times, calculations shows that Z is much smaller than 1.0, 
especially around the Dirac point. This unusual feature directly reflects that a Bloch states 
around the Dirac point couples strongly with other nearby Bloch states. Therefore, the Bloch 
states would not be a good starting point to perform perturbative calculations. In fact, the 
similar feature has been observed in semiconductor alloys and has been termed non-Bloch 
nature of alloy states 30(]. This is why the SCBA method cannot produce accurate enough 
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FIG. 3: (Color online) (a) Effective energy dispersion relation for ideal graphene and disordered 
graphenes for riimp/N = 10% and V = 0.5 (filled squares) and V = 2.0 (filled circles), (b) The 
ratio Vg/vp (filled symbols) and Vp/vp (open symbols) as a function of wavenumber k (in units of 
1/a). 

results. This nontrivial feature could lead to intriguing transport and optical properties. For 
example, linear dispersion relation leads to En = ±VpV2ehnB {n = 0, 1,2...) LL series. 
Since the LLs come from the quench of electron's kinetic energies, lower group velocity 
near the Dirac point means the energy gap between n = LL and n = ±1 LL should be 
substantially reduced by the disorders while that between n = ±1 LLs and n = ±2 LLs will 
be affected less by disorder. As a result, one should expect the ratio (£'2 — E^i)/{Ei —En) is 
substantially larger than v^+ 1- This is indeed what was observed in a recent experiment 311 ] 
where {E2 — E^i) / (Ei — Eq) is found to be greater than -\/2 + l, by 5 ~ 10%. This observation 
implies that their samples are quite disordered. 

Conductivity. — Assuming that the Einstein relation for the Dirac electrons is valid, one 
can obtain the conductivity from a = e^Dp, where e is the electron charge, D = v'lTtr/2 
is the diffusion coefficient, and Ttr denotes the transport relaxation time and satisfies the 
relationship Ttr = 2re for the point-like disorder. [l7| Fig. 4 is the conductivity calculated 
from transport relaxation time Ttr = 2re in Fig. 2, the group velocity Vg in Fig. 3, and 
the density of states obtained by p{E) = —^Jq" -J^K^A^k, E) as shown in inset to Fig. 
4. For very weak disorder, the conductivity is sub-linear in the Fermi energy. Increasing 
the disorder strength, the conductivity weakly depends on the Fermi energy or the carrier 
density, in agreement with the prediction of Boltzmann theory for short-range scatterers 
[l7|. At the Dirac point, by combining the self-energy function with Kubo formula, the 
conductivity is predicted to have the universal value crxx{0) = 4e^/7r/;,|i32|]. However, our 
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FIG. 4: (Color online) Conductivity as a function of charge density. The model parameters are 
nimp/]^ = 10% and Vq =0.5, 1.0 and 2.0. Inset: Density of states for the same parameters. 

calculations find that ^^^^^(O) takes the values 8.2e^//i, 5.9e^/h and 2.7e^//i for a = 0.0046, 
0.018 and 0.07, respectively, and shows non-universal behavior. 

Summary. — In conclusion, we studied point-like disorder effects on the one-electron 
properties of graphene. The exact ensemble-averaged Green's function is obtained from a 
large-scale real-space calculation. Through the analysis of self-energy and spectral functions, 
we conclude that the single-particle lifetime reduction and the linear dispersion relation are 
modified by the hybridization of the Bloch states. Furthermore, we studied the diffusion 
transport properties by using our exact self-energy and the Einstein relation. Our approach 
is very general and robust, thus is apphcable to many other disordered systems. 
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